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in : A microscopic theory of a non-dissipative drag in a two-component superfluid Bose gas is devel- 

■ oped. The expression for the drag current in the system with the components of different atomic 

' masses, densities and scattering lengths is derived. It is shown that the drag current is proportional 

to the square root of the gas parameter. The temperature dependence of the drag current is studied 
(~| . and it is shown that at temperature of order or smaUer than the interaction energy the tempera- 

'~{ ' ture reduction of the drag current is rather smalL A possible way of measuring the drag factor is 

proposed. A toroidal system with the drag component confined in two half-ring wells separated by 
. two Josephson barriers is considered. Under certain condition such a system can be treated as a 

Bose-Einstein counterpart of the Josephson charge qubit in an external magnetic field. It is shown 
that the measurement of the difference of number of atoms in two wells under a controlled evolution 
of the state of the qubit allows to determine the drag factor. 



I 



PACS numbers: 03.75.Kk, 03.75.Lm, 03.67.Lx 



INTRODUCTION 



'"^ Macroscopic quantum coherence manifests itself in many specific phenomena. One of them is a non-dissipative drag 

' that takes place in superfluids and superconductors. The non-dissipative drag, also known as the Andreev-Bashkin 
O ' effect, was considered, for the first time, in Ref. where a three velocity hydrodynamic model for '^He-'*He superfluid 
mixtures was developed. It was shown that superfluid behavior of such systems can be described under accounting 
the "drag" term in the free energy. This term is proportional to the scalar product of the superfluid velocities of 
two superfluid components. A similar situation may take place in mixtures of superfluids of 5z = -f 1 and 5*2 = —1 
^ ' pairs in liquid '^He in the A-phase Among other objects, where the non-dissipative drag may be important, are 
, neutron stars, where the mixture of neutron and proton Cooper pair Bose condensates is believed to realize jj,!^. 
' The possibility of realization of the non-dissipative drag in s^erconductors was considered in . The non-dissipative 
" drag in bilayer Bose systems was treated microscopically in [q, for a special case of two equivalent layers of charged 
■ bosons. The case of a bilayer system of neutral bosons was studied in Q in the limit of small interlayer interaction. 
' The most promising systems where the non-dissipative drag can be observed experimentally are two-component 
, alkali metal vapors. In such systems the interaction between atoms of different species is of the same order as the 
' interaction between atoms of the same specie and the effect is expected to be larger than in bilayers. In Bose mixtures 
the components are characterized by different densities, different masses of atoms and different interaction parameters. 
In this paper we consider such a general case and obtain an analytical expression for the drag current for zero and 
finite temperatures. 

In the system under consideration the drag force influences the dynamics of atoms in the drag component in the 
Q ■ same manner as the vector potential of electromagnetic field influences the dynamics of electrons in superconductors. 
^ In particular, in neutral superfluids with Josephson links the drag effect may induce the gradient of the phase of the 
^ ' order parameter in the bulk and, as a consequence, control the phase difference between weakly coupled parts of the 
. ^ , system. Therefore, on can expect that the effect reveals itself in a modification of Josephson oscillations between 
' weakly coupled Bose gases. In this paper we discuss possible ways for the observation such a modification. We 
, consider the Bose gas confined in a toroidal trap with two Josephson links. In the Fock regime Q the low energy 
dynamics of the system can be described by the qubit model of general form (the model, where all three components 
of the pseudomagnetic field can be controlled independently). The parameters of the qubit Hamiltonian depend on 
the drag factor. The measurement of the state of the qubit under controlled evolution allows to observe the effect 
caused by the non-dissipative drag and determine the drag factor. In this paper we consider two particular schemes 
of the measurement. In the first scheme one should determine the time required to transform a reproducible initial 
state to a given final state. In the second scheme the geometrical (Berry) phase should be detected. 

In Sec. mthe microscopic theory of the non-dissipative drag in two-component Bose gases is developed. In Sec. IIIII 
a model of the Bose-Einstein qubit subjected by the drag force is formulated and the schemes of measurement of the 
drag factor are proposed. Conclusions are given in Sec. IIVI 
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II. NONDISSIPATIVE DRAG IN A TWO-COMPONENT BOSE SYSTEM. MICROSCOPIC 

DERIVATION 

Let us consider a uniform two-component atomic Bose gas in a Bose-Einstein condensed state. We will study the 
most general situation where the densities of atoms in each component are different from one another (ni n2), the 
atoms of each components have different masses {mi ^ 1112) and the interaction between atoms is described by three 
different scattering lengths (an 7^ 022 7^ 012) • The Hamiltonian of the system can be presented in the form 
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i,i'=l,2 



where 



E:= J dV^[V*+(r)]V*.(r) (2) 

is the kinetic energy, 

E^'^ I d\4>+{r)^+{r)Ju'■^Ar)■^^{r) (3) 

is the energy of interaction, = Airh^ an / nii and 712 = 2iTh^{mi + m2)ai2 / {^1111712) are the interaction parameters, 
and fii are the chemical potentials. 

For the further analysis it is convenient to use the density and phase operator approach (see, for instance, [Tol[ll| ). 
The approach is based on the following representation for the Bose field operators 



^'i(r) = exp [i(pi{v) + iipi{v)] \Jni+ ni{r), (4) 



^'+(r) = y^ni + ni{r) exp [-i(pi{r) - i(pi{r)] . 



(5) 



where rii and tpi are the density and phase fluctuation operators, ipiir) are the c-number terms of the phase operators, 
which are connected with the superfluid velocities by the relation — hVipi/rrii. In what follows we specify the case 
of the superfluid velocities independent of r. 

Substituting Eqs. I@J, (O into Eq. JQ) and expanding it in series in powers of fii and \7(pi we present the Hamiltonian 
of the system in the following form 



H — Hq + H2 



(6) 



In ® the term 



, 1=1,2 



1 2 O'ii 2 



- 7i2nin2 



(7) 



does not contain the operator part. Here V is the volume of the system. The minimization conditions for the 
Hamiltonian Hq yield the equations 



1 



^1^1 + lii^i + 7i2n-3-i - Mi = (i = 1, 2). 



(8) 



Under the conditions (jSJ the terms, linear in the density fluctuation operators, vanish in the Hamiltonian. Taking 
into account the V(niV(y3i(r)) = 0, we find that the terms, linear in the phase fluctuation operators, vanish in the 
Hamiltonian as well. 

The part of the Hamiltonian quadratic in V(^i and hi operators reads as 



2m. 



(Vn.(r))^ 

4n,- 



+ ^(n,(r)V<p^(r) + [V^,(r)]n,(r) 
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|^([Vn,(r)]V^.(r) - [V(^.(r)]Vn,(r)) + ^(n,(r))' I + 7i2ni(r)n2(r) 



(9) 



The quadratic part of the Hamiltonian determines the spectra of the elementary excitations. Hereafter we will 
neglect the higher order terms in the Hamiltonian ©. These terms describe the scattering of the quasiparticles and 
they can be omitted if the temperature is much smaller than the temperature of Bose-Einstein condensation. 

Let us rewrite the quadratic part of the Hamiltonian in terms of the operators of creation and annihilation of the 
elementary excitations. As the first step, we use the substitution 



*^ , V ^ik 

k 



(10) 



kr 



[6,(k)-6+(-k)] 



(11) 



where operators bf, bi satisfy the Bose commutation relations. Here etk = fi^k^ /2mi is the spectrum of free atoms, 



and 



Elk = \J fik{eik + "^luni) 



(12) 



is the spectrum of the elementary excitations at 712 — and v.^ = 0. The substitution (|10|l . reduces the 
Hamiltonian to the form quadratic in hf and hi operators: 



£,(k)(6+(k)6,(k) + i') -ie,. 



ik 

+ Y.3^ [''^(k)^2(k) + 6i(k)62(-k) + h 



(13) 



Here 



and 



fi(k) = Eik + hkvi 



gk = 712 



eifc£2fc?^i>t2 
EikE2k 



(14) 
(15) 



The Hamiltonian (|13|l contains non-diagonal in Bose creation and annihilation operator terms and it can be diago- 
nalized using the standard procedure of u-v transformation ,13.] . The result is 



^2 = E E ^^(k) ( /3+(k)A(k) + 1 ) - i ^ e.. , 



where (3^ {k) and f3\{k) are the operators of creation and annihilation of elementary excitations. 
The energies ^^A(k) satisfy the equation 



det 



where 



'fi(k) 

Qk 




B 



-El 




^ ) = 
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+ E1) 







gk 


£i{- 


k) 










f2(k) 


gk 




£ 







gk\ 








gk 





gk 





\gk 





/ 



gk 




(16) 



(17) 



(18) 



(19) 
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and I is the identity matrix. 

The densities of superfluid currents in two components can be obtained from the relation 



. _ 1 

~ y 9^' 



(20) 



where F is the free energy of the system. Here the quantity is defined as the density of the mass current. 
The free energy of the system, described by the Hamiltonian |(BJ, is given by the formula 



J2 ^A(k)- J2 ^''^ 

\=a,0 1=1,2 



E 1- 



1 — exp 



gA(k) \ 



(21) 



The second term in H21|l is the energy of the zero-point fluctuations and the third term is the standard temperature 
dependent part of the free energy for the gas of noninteracting elementary excitations. 

We specify the case of small superfluid velocities (much smaller than the critical ones). In this case the currents 
can be approximated by the expressions linear in Vi . To obtain these expressions we will find the free energy as series 
in Vi, neglecting the terms higher than quadratic. 

At vi = V2 = the equation H17|) is easily solved and the spectra are found to be 



E. 



a(l3)k 



2k 



± 



1/2 



47?2"i"-2eifee2fe 



(22) 



As required in the procedure [I^, we take positive valued solutions of Eq. H17|) . The energies (|22|l should be real 
valued quantities. This requirement yields the common condition for the stability of the two-component system: 
7^2 ^ 7ii722- If this condition were not fulfilled, spatial separation of two components (at positive 712) or a collapse 
(at negative 712) would take place. 

At nonzero superfluid velocities we present the solutions of Eq. p7|l as series in v^: 



i;ikv2 ( 1 - 



27i2'^i"-2eifee2fe 

Eak {eIj^ - El^ 



(kvi — kv2) 



(23) 



Spik) = Epk + 2^kvi I 1 



^ife 



E2k 



p2 _ p2 
^ak ^pk 



-ftkv2 1 



Elk - E2k 



^Ik 



^Ik 



2-il2nin2eiki2k(El^ + 2>El^ 

^ —^h^ (kvi - kv2)^ 



Epk 



{^Ik 



Epk 



(24) 



Note that at vi = V2 = v the spectra (|23|) . (|24() are reduced to common expressions for the energies of quasiparticles 
in a moving condensate: f a( /3)(k ) — Ea(i3)k + ^kv. 

Using Eqs. H21|l . (|23|l and H24|) we obtain the following expression for the free energy 



F = Fo 



V 



(pi - Pnl)Vi + (p2 - Pn2)v2 - Pdr(vi - V2)^ 



(25) 



where Fq does not depend on Vj. In H25(l pi = rriini are the mass densities, the quantities 



Pnl 



mi 



E^ife 

k 



IK 



dEpk El^~Ej^ \dEak dEpk 



(26) 



Pn2 



1712 



E^2fc 



dN, 



ak 



dN, 



pk 



Elk 



Elk ( dNak 



dN, 



Pk 



dEak 



dE^ 



pk 



Et 



Epk \d-Eak dEpk 



(27) 
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describe the thermal reduction of the superfiuid densities, and the quantity 



Pdr 



-V^2}^ 



1 + jVgfc + Ni3k 

{Eak + El 



Nr. 



2EnkE, 



akJ^pk 



dEak 



dNpk 
dE. 



■(3k 



(28) 



which we call the "drag density," yields the value of redistribution of the superfiuid densities between the components. 



In Eqs. S-lEHIl Nr^k 
Using Eqs. 



[exp(ii'Q(^)i./r) — 1] ^ is the Bose distribution function, 
we arrive to the following expressions for the supercurrents 



jl = (Pl - Pnl - Pdr )vi + pdrV2, 



(29) 



j2 = [p2 - Pn2 - Pdr)^2 + PdrVl . 



(30) 



One can see that at nonzero pdr the current of one component contains the term proportional to the superfiuid velocity 
of the other component. It means that there is a transfer of motion between the components. In particular, at wi = 
the current in the component 1 (ji = Pdr^2) is purely the drag current. Since pdr is the function of 7^2 (see Eqs. (|28|) 
and 12211 ) the drag current does not depend on the sign of the interaction between the components. 

Eq. (|28|l is the main result of the paper. This equation yields the value of the drag for the general case of 
two-component Bose system with components of different densities, different masses of atoms, different interaction 
parameters, and for zero as well as for nonzero temperatures. Moreover, this equation is valid not only for the point 
interaction between the atoms, but for any central force interaction. In the latter case the interaction parameters 
7ifc in Eq. H28|) and in the spectra 1)22(1 . (|12|l should be replaced with the Fourier components of the corresponding 
interaction potentials. 

To estimate the absolute value of the drag we, for simplicity, specify the case mi = m2 = m, that is realized when 
two components are two hyperfine states of the same atoms. 
At T = Eq. 1211) is reduced to 



4m , 
Pdr^— I de 



7?2'^i^2i^(e)e^/^ 



3 Jo ^(e + wi){e + W2) (\/e + wi + a/c + ^2)" 



(31) 



where 



is the density of states for free atoms, and 



i'(e) ~ —pz -v/e 



wi{2) = 7ii"i + 722"-2 ± \J iliini - 722*^2)^ + 47^2"1"2- 



The integral in ((31(1 can be evaluated analytically. To present the answer in a compact form it is convenient to 
introduce the dimensionless parameters 



*12 



011022 



and 



n2a22 



^2022 
niaii 



(0 < ?7 < 1 and K > 2) 

Using these notations we have 



Pdr 



VPiP2\/ nial^n2al2^F{K, 77), 



(32) 



where 



256 



45V27r 



(k + 3^1 - 



^ K.+ yf^- 4 + 4ry + ^k - ^/k?' - A + Arj 



(33) 
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FIG. 1: Dependence of the "drag density" on the temperature. 



Direct evaluation of Eq. H33(l shows that at allowed rj and k (0 < 77 < 1 and k > 2) the factor F{k, rj) is almost the 
constant (the range of variation of F is [0.7 ~ 0.8]) and one can neglect the dependence of F on the parameter of the 
system. 

At oiini = 022772 we obtain from 1)32(1 the following approximate relation 



Pdr 



1 



'■12 



aiia22 



1 



2 011022 



12 'n242. 



(34) 



If the density of one component is much larger than of the other and an ^ 022, the "drag density" is approximated 
as 



Prfr ~ 0.8^1-2-^/ 772022 at 771^ 772, 



Pdr ~ 0.8p2-p- V '^lofi at 772 < 771. 



(35) 



One can see that the "drag density" is proportional to the square root of the gas parameter. It means that the drag 
effect is larger in "less ideal" Bose gases. 

The temperature dependence of the " drag density" at small T can be evaluated analytically from Eq. (|28() using 
the linear approximation for the spectra of the excitations. It yields pdr{T) — pdr(0)(l — aTT^/T^), where Tq = 
^7ii77i722772 and the factor ut is positive. Numerical evaluation of the sum over k in Eq. ((28|l shows that the 
analytical approximation is valid only at T <C Tq. At T > Tq the "drag density" decreases much slower under increase 
of the temperature. As an example, the dependence of pdr{T) at 77i = 772 = 77, 711 = 722 = 7 and rj = 0.5 is shown in 
Fig. CI 

Now let us discuss how the drag effect can reveal itself in a real physical situation. If one deals with the stationary 
superflow one implies that it is the circulating superflow, e.g., the tangential superflow in a hole cylinder. In such a 
case the superfluid velocities satisfy the Onsager-Feymnan quantization condition 



2TrnNj 



(36) 



where the vorticity parameters Ni are integer. Then, the drag effect can be understood as the appearance of the 
circulating current in the drive component (e.g. specie 1), when the circulation of the superfluid velocity of the drive 
component (e.g. specie 2) is fixed (A^2 = const). The current of the specie 1 1(29(1 depends on the superfluid velocities 
of the both species and if the superfluid velocity of the drag component directed anti-parallel to the superfluid velocity 
of the drive component the current of the drag component might vanish. But since the velocities are quantized it may 
happen only under certain special conditions (see below) . The superfluid velocity of the drag component is determined 
by that at fixed A'^2 the free energy 1(25(1 has a minimum with respect to discrete values of Vi = HNi/ {miR) (where the R 
is the radius of the contour in ((36() ). Depending on the value of the parameter a = {pdr/ {P2 — Pn2 — Pdr)){'mi/m2)N2 
several possibilities can be realized. At |a| < 1/2 the minimum of the energy 1(25(1 corresponds to A^i = (and 
Vi =0). In this case the current of the drag component is directed along the drive current and it is proportional 
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to the drag density. At |a| = p (p is natural) the value A^i = ~p minimizes the energy. In this case two terms in 
Eq. (|29|) compensate each other and the current in the drag component vanishes. At half-integer a the degenerate 
situation takes place: two state (with co-directed currents, and counter-directed currents) have the same energy. At 
1/2 + p < \a\ < p the state with counter-directed currents gains the energy and at p < \a\ < p + 1/2 the co-dirccted 
currents are energetically preferable. In the latter two cases the nonzero vorticity of the drag component (A^i ^ 0) is 
also induced. This behavior is analogous the behavior of a superconducting ring in a magnetic field. We note that 
since pdr ^ P2, the most realistic case is \a\ < 1/2 when the simple picture of the transfer of part of the motion from 
the drive to the drag component takes place. 

In this study we have concentrated on the analytical derivation of the drag effect in the uniform Bose gases. The 
consideration of the non-uniform case requires the solution of the eigenvalue problem for the elementary excitations 
in the two-component Bose gas in the external potential. But even for the simplest case of a spherically symmetric 
trap this problem can be solved analytically only in the long-wavelength limit and the Thomas-Fermi approximation 
|l2| (the spectrum of elementary excitations in one-component Bose gases was obtained analytically for a number of 
potentials but also in the same Hmit 0,0,E,0|)- Since the main contribution to the drag density comes from 
the excitations with the wave vectors of order of the healing length (see (jSHl), the rigorous analysis of the drag effect 
can be done only numerically. Nevertheless, in the Tomas-Fermi situation the drag effect can be evaluated basing on 
the following arguments. When the linear size of the Bose cloud is much larger than the healing length, the spectrum 
of the excitations at the wave vectors of order or higher than the inverse healing length is well described by the quasi 
uniform approximation. Therefore, the drag effect can be described by the same equations, as in the uniform case 
with the only modification that the quantities ni and 712, and, correspondingly, pi, pni, Pdr and ji in Eqs. I|26|I - I|3U|) 
are understood as functions of coordinates. 

At an arbitrary symmetry of the trap potential the superfluid velocity of the drag component cannot be equal to 
zero in each point. Indeed, in general case of space dependent pi, pni, and pdr the velocity field V2(r) cannot satisfy 
two independent continuity conditions V[(p2 — Pn2 — Pdr)^2] = and V{pdr'^2)- To analyze this case one should find 
the velocity fields vi(r) and V2(r) that satisfy the continuity conditions and the quantization conditions. To illustrate 
this point let us consider a simple example of a trap having the shape of a hollow cylinder with the densities that 
depend only on the polar angle <j). We will seek the velocity fields that do not have radial components. Then, the Eqs. 
(|29|l . H30() . written for the tangential components of the currents and the velocities, can be presented in the matrix 
form 

where 



^ f PslW - Pdr{(f>) Pdri(t)) \ /ggN 

V P<ir(0) Ps2{(t)) - Pdr{4>) ) 

with Psi{(t>) — Pi(0) — Pni{<P)- Due to the continuity conditions the current ji and j2 in (|37|1 do not depend on <j). 
According to Eq. (|37|) the velocities vi (r, <j)) and V2 {r, cj)) are connected with the currents by the equation 

^i("'<^)^-i?-i('^>!"n (39) 



V2{rA) J \h{r) 

Integrating Eq. (|39|l over (p and taking into account the quantization conditions H3t)|l we obtain the equation for the 
currents 

f(n{r)\_2nh^N^/m^\ 

^ \32{r))" r yN^/m^)^ ^^"^ 

where 

/ f dS Ps-i-Pd. _ f^^^ d6 2£fc \ 

rf _ ( Jo ^ PBlPB2-Pdr-(Psl+PB2) JO ^ P si P b2 - Pdr {P si + P s2} \ /At\ 

- f'" deb f'" dd P-^-^"- j ^ ' 

\ Jo ^ PslPs2-Pdr(Psl+Ps2) Jo ^ PslPs2-PdripBl+Ps2) / 

If a given vorticity of the drive component A'2 is not very large the minimum of energy is reached at A^i = 0. In the 
latter case the solution of Eq. H4()|l in the leading order in pdr yields the following expression for the current of the 
drag component 

r27V ,1 pdr{4>) 

j ir) ~ 27r^^2 Jo ^'Pp.iWp.sW .^2) 
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FIG. 2: Schematic shapes of Bose clouds for the drag (top figure) and drive (bottom figure) components. The drag component 
is situated inside and overlaps with the drive component. 



One can see that if at some 4> the density psi has a sharp minimum the first factor in denominator in Eq. (|42|l becomes 
large. On the other hand, the integral in numerator is not very sensitive to lowering of psi (see Eqs. (jSHJ)- Thus, in a 
system with a "bottle neck" in the drag component the drag current decreases strongly and the main consequence of 
the drag effect is the emergence of the gradient of the phase of the order parameter of the drag component. Similar 
situation takes place in a system with a weak link. The latter case is analyzed in the next section. In the uniform 
case Eq. (|^ is reduced to ji = PdrV2- 

To complete the discussion we emphasize that the crossed term hr the free energy (|25|l (and, consequently, 

the drag terms in the currents H29|) . H3U|l ) comes only from the second and third terms in Eq. H21|) . Consequently, the 
drag effect considered in this paper is solely by the excitations. At the mean field level of approximation (which can 
be also formulated in terms of the Gross-Pitaevsky equation) the effect does not appear, while the coupling between 
the components is also present at that level of approximation. We would note that at the mean field level the drag 
effect of another type may emerge. That effect takes place in the case when one of the species is subjected by an 
asymmetric rotating external potential (see, for instance, 17], where such an effect has been studied with reference 
to the system of two coupled traps). 



III. MODEL OF BOSE-EINSTEIN QUBIT WITH EXTERNAL DRAG FORCE 

It is known that Bose systems in the Bose-Einstein condensed state may demonstrate Josephson phenomenon j^. 
It this paper we consider the external Josephson effect that takes place in two-well Bose systems. It was shown in 
that in such systems one can realize the situation, when two states, that differ in the expectation value of the 
relative number operator, can be used as qubit states. 

To include the drag force into the play we consider the following geometry. Let our two-component system is confined 
in a toroidal trap and the Bose cloud of the component 1 (the drag component) is situated inside and overlaps with 
the Bose clouds of the component 2 (the drive component). Such a situation can be realized if I712I < min(7ii,722). 

Deforming the confining potential one can cut the drag component into two clouds of a half-torus shape (separated 
by two Josephson links) leaving the Bose cloud of the drive component uncutted ( Fig. EJ. In what follows we use 
the following notations: Rt is the large radius of the toroidal trap, rti and rt2 are the small radiuses of the toroidal 
Bose clouds of the drag and the drive components, correspondingly. 

Rotating this trap one can excite a superflow in the drive component. After the rotation be switched off there 
will be a circulating superflow in the drive component and no superflow in the drag component (at negligible small 
Josephson coupling). The superfluid velocity of the drive component is 

«2 = 43 

In (|43|l . we imply that Rt 3> rt\, rti and neglect, for simplicity, the effect caused by a dependence of rt2 on the polar 
angle. 
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Since ji = 0, the phase gradient Vipi should be nonzero to compensate the drag effect, 
the (f) component of the phase gradient is given by the relation 



iV2 

Rf. 



fdr = -fdr{^V2)<t,, 



In the polar coordinates 



(44) 



where 



dr 



Pdr 



"^2 Psl - Pdr 



(45) 



The quantity fdr yields the ratio between the phase gradients in the drag and the drive components in the situation 
when the drag component is in the open circuit (i.e. the current cannot flow in the circuit). We call this quantity the 
drag factor. 

We imply that rn and rt2 are much larger than the healing lengths that allows to describe the drag effect in quasi- 
uniform approximation. For definiteness, we specify the case of pi <C P2 and p2 ~ const in the overlapping region. In 
this case one can neglect the space dependence the drag factor (see Eqs. (j^ ). 

At nonzero Josephson coupling the current ji can be nonzero, but it cannot exceed the maximum Josephson current 
jm- Relation H44|l remains approximately correct at nonzero Josephson coupling, if an inequality jm <C hpi/ (jniRt) 
is satisfied. Here we specify just such a case. It is important to emphasize that we consider the situation, when there 
is only the external Josephson effect between two half-torus traps, and there is no internal Josephson effect between 
the drag and the drive species. 

The drag force can be considered as an effective vector potential Adr = — ft/drV(^2 (in units of e = c = 1) that 
corresponds to an effective magnetic flux ^dr = — 27rfi./dr-/V2- Thus, our Bose system is similar to the Cooper pair box 
system that implements the Josephson charge qubit with the Josephson coupling controlled by an external magnetic 
flux Tg]. To extend this analogy we formulate the model of the Bose-Einstein qubit subjected by the drag force. In 
what follows we use the approach of Ref . . 

In the two mode approximation the Bose field operators for the drag component can be presented in the form: 



4'i(r,i) 



E 

l=L,R 



ai{t)^i{r-ri), 4'^(r,t) 



l=L,R 



(i)*r(r-rO, 



(46) 



where 



and a^ji) are the operators of creation and annihilation of bosons in the condensates confined in the 



left (right) half-torus, and l, ^! b. are two almost orthogonal local mode functions 



/, /' — R 



that describe the condensate in the left and right traps '2(?| . 

Substituting H46|l into Hamiltonian we obtain the following expression for the parts of the Hamiltonian that 
depends on the operators and a/: 



(47) 



l=L,R 



with 



Ki = / dV*;- 



2m 



(48) 



7ii 
2 



(49) 



J ^ d 



Zm 



(50) 



The functions "^l and contain the phase factors e*'''^*^'"^ and e*'^^'-'"-', where the phases satisfy Eq. (|44ll . Taking 
these factors into account, one can choose the following basis for the one mode functions 



*L(fl)(r) = |*L(fl)(r)|exp [-iN2fdr(l>UR)(r)] , 



(51) 
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where 4>L,<t>R are the polar angles counted from the centers of L and R half-torus, correspondingly (see Fig. The 
angles (f>L(R){i^), defined as shown in Fig. 13 satisfy the relation 

0fl(rA) - (I^LirA) = (f>L{rB) - <l>R{rB) = tt, (52) 

where ta and are the radius- vectors of Josephson links. 

Substituting (|51(l into Eq. H50(l . using Eq. Il52|l and taking into account that the functions ^'l and \E'/j overlap in 
a small vicinity of A and B links, we obtain the following expression for the Josephson coupling parameter: 



J (Ja + Jb) cos (^-^^ + i{JA - Jb) sin ' 



where $o — 27r/i is the "flux quantum" and 



Ja(b) ~ / d' 



j3 " 



r 



Zm 



(53) 



(54) 



Here Va and Vb are the areas of overlapping of two one mode functions at links A and B, correspondingly. 
Considering the Hilbert space in which the total number operator 

N ^ aj^QL + a^aji (55) 

is a conservative quantity {N = N) we present the Hamiltonian (|47|l in the following form 

Ha = EcifiRL - rigf + {Jaj^aii + h.c.) + const, (56) 

where 



UR. = <^---t^^ (57) 



is the number difference operator, 

Ec = Xr + Xl (58) 

is the interaction energy, and the quantity 

ng^-^[KL-KB^ + {N-l)iXL-\R)] (59) 
describes an asymmetry of L and R half-tore. 
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In what follows we imply that the system is in the Fock regime Q {\J\N <^ Ec) and use the number representation 

\nRL) = \nB.,nL) = ly + '^RL, y ~ "-fli)- 

In this representation the first term in (|56|l is diagonal. The second term in H56(l can be considered as small nondiagonal 
correction. But if Ug is biased near one of degeneracy points 

for even iV (gp) 
^ \m for odd ^ ^ 

(where M is integer and \M\ < N/2), the second term in (|56|l results in a strong mixing of two lowest states 
(I T) — \ndeg + 1/2) and | \) = \ndeg — 1/2)) and the low energy dynamics of the system can be described by a 
pseudospin Hamiltonian 



where di are the Pauli operators, and 



Heff = ^O-^ - -^CTy - ^Cr^, (61) 



fix = -(Ja + Jb)\I{N + 1)2 - 4n2^g cos ( vr^ 



= -{J A - Jb)\I{N + 1)2 - 4^2^^ sin ( TT^ 



flz = 2Ec{ng - Udeg) (62) 



are the components of the pseudomagnetic field. In experiments one can control the parameters n^, J a and Jb 
independently and, consequently, the pseudomagnetic field fl{t) can be switched arbitrary. It mean that Eq. 
represents the standard Hamiltonian of the qubit system. The parameters of the qubit (|61|l depend on the "drag 
flux" Therefore, one can determine its value from the measurement of the state of the system after a controlled 
evolution of a certain reproducible initial state. 

Let us consider two possibilities. For definiteness, we specify the case of odd N and the degeneracy point Udeg = 0. 

If the Josephson coupling are switched off and rig is switched on to some positive value (much less than unity) the 
system is relaxed to the state IV'm) = I ?)• This state can be used as the reproducible initial state. The quantity 
should be measured is the expectation value of the number difference operator. In the initial state the expectation 
value of this operator is url = 1/2 

When the system is switched suddenly to the degeneracy point Ug = and the Josephson couplings are switched 
on for some time r the initial state evolves to another state with another url- 

If one sets Ja = Jb = J the result of evolution {\ipf) = U\^pin)) is described by the unitary operator 

jj ,\ _ f cos(aiT) — isin(aiT) 
I — isin(Q!ir) cos(Q!ir) 

where ai = {J/h){N + 1) cos(7r$dr/5'o)- One can see that at time of evolition r = ti = 7r/(4|Q:i|) the expectation 
value of the number difference operation will be equal to zero. 

For the case Ja — J and Jb — the operator of evolution reads as 



U2{r) 



cos(a2T) —ie ■"r*<j>-/<i'o sin(a2T) 

_^g27r$dr/*o sin(a2T) cos(a2T) 



with 02 = {J/2h){N + 1). Respectively, the expectation value rifjL will be equal to zero at r = r2 = 7r/(4a2) 

The ratio T2/T1 — \ cos(7r$dr/*I'o)|/2 depends only on ^dr and the quantity ^dr can be extracted from the measure- 
ments of Ti and T2 . It is important to note that to provide this scheme one should control only the ratio of J a and 
Jb, but not their absolute values. 

Another possibility can be based on detection of the Berry phase Eq. H61|) contains all three components of 
the field Jl and they can be controlled independently. The general scheme of detection of the Berry phase in such 
a situation was proposed [S^. A concrete realization of this scheme in the Josephson charge qubit was described in 
[2^ . Here we extend the ideas of jU|2^ to the case of the "dragged" Bose-Einstein qubit. 

We start from the same initial state and switch to Ja = Jb ~ J and ng = 0. The initial state | t) can be presented 
as the superposition of two instantaneous eigenstates \ea) = (| T) + I i))/V2 and \eb) = (| T) — I i))/^' 

= ^(|e.) + \et)). (63) 



12 



An adiabatic cyclic evolution of the parameters of the Hamiltonian (|61(l results in appearance of the Berry phase in 
the \ea) and |ef,) eigenstates, if the vector Q, subtends a nonzero solid angle at the origin. 

Let us consider the following 4 stage cyclic adiabatic evolution starting from the point Ja = Jb = J and Ug — 0: 
1 - Js is switched off; 2 - J^i is switched off and simultaneously Ug is switched to Ugi > 0; 3 - rig is returned to the 
same degeneracy point {ug ~ 0) and Jb is switched to Js = J; 4 - Ja is switched to Ja = J (all switches should be 
done slowly: h\d^l/dt\ < il^). 

After such an evolution the system arrives at the state 

|^™) = -^(e^*"+'-'|e,)+e^^-^^|e,)), (64) 

where 7 = n^dr/^o is the Berry phase (equals to half of the solid angle subtended by fJ) and 5a, Sb are the dynamical 
phases. 

Elimination of the dynamical phases can performed by swapping the eigenstates (7r-transformation) and repeating 
the same cycle of evolution in a reverse direction (see |2a|'). 

The TT-transformation can be done by fast switching off the Josephson coupling and switching on Ug — ng2 > 
during the time interval t^^ = tm / {2Ecng2). After the 7r-transformation the state becomes 

= --^(e'*"+*^|e,) + e'^^-'^le,)). (65) 

After the cyclic evolution in the reverse direction we arrive to the state 

\^f) - --J=e^(^"+^^)(e^^^|e,) + e-^^^je,)). (66) 

One can see that the expectation value of the number difference operator in the final state H66|l uhl — cos(47)/2 — 
cos(47r$dT-/^o)/2 depends only on $dr and the measurement of this difference allows the determine the value of the 
"drag flux". 

Thus, the measurements of relative number of atoms in left and right condensates under controlled evolution of the 
state of the system allows to observe the non-dissipative drag and determine the drag factor (if the vorticity of the 
drive component is known). 

IV. CONCLUSIONS 

We have investigated the non-dissipative drag effect in three-dimensional weakly interacting two-component super- 
fluid Bose gases. The expression for the drag current is derived microscopically for the general case of two species 
of different densities, different masses and different interaction parameters. It is shown that the drag current is pro- 
portional to the square root of the gas parameter. The drag effect is maximal at zero temperatures and it decreases 
when the temperature increases, but at temperatures of order of the interaction energy the drag current remains of 
the same order as at zero temperature. 

We have considered the toroidal double-well geometry, where the non-dissipative drag influences signiflcantly on 
the Josephson coupling between the wells. In the system considered the drag force can be interpreted as an effective 
vector potential applied to the drag component. The effective vector potential is equal to A^^ = —hfdr^fdrv (in units 
of e = c = 1), where ipdrv is the phase of the drive component, and fdr is the drag factor. In the toroidal geometry 
the effective vector potential can be associated with an effective flux of external field ^dr = '^T^fifdrN^, where iV„ is 
the vorticity of the drive component. In the Fock regime the system can be considered as a Bose-Einstein counterpart 
of the Josephson charge qubit in an external magnetic field. The measurement of the state of such a qubit allows to 
observe the drag effect and determine the drag factor. 
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